Nonlinear optical responses and gap-soliton lattices in layered dielectric materials. 
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A new approach to the nonlinear reponse of a multiplayer structure consisting of alternating di- 
electric materials with intensity-dependent dielectric constants is considered in both analytical and 
numerical viewpoints. The nonlinear feedback of the structure gives rise to an artificial potential 
acting as a mean-field feedback on each single dielectric material due to other layers. For superlat- 
tices resulting from a regular periodic arrangement of an infinite number of single dielectric layers, it 
is shown that an analytical expression of this mean-field feedback can be obtained and waveshapes 
of the associate gap solitons are derived. The case of a multilayer with dispersive interlayer inter- 
actions is considered numerically and leads to strongly localized gap solitons. Implicatons of the 
artificial soliton structure as for possible theoretical modelling of soliton-compression phenomena 
are discussed. 

PACS numbers: 42.55.Tv, 42.65.Sf, 42.65.Tg 
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Superlattices with Bragg reflectors are generally not 
accessible to any light frequency. In multilayer struc- 
tures consisting of a periodic arrangement of single di- 
electric materials, the dependence of dielectric constants 
on the light intensity affects the intimate structure of 
materials leading to the formation of defects which are 
rcponsiblc for these apparent gap excitations. Chen and 
Mills pointed out a mechanism so-called stop-gap ra- 
diation transmission following which the light could still 
propagate within the superlattice due to partial closing of 
the gap defect. Later, it was suggested^, that a suit- 
able arrangement of single dielectric materials composing 
a superlattice structure too may favor gap-soliton propa- 
gation. In recent years, a great deal of interest has been 
devoted to understanding nonlinear features of photonic 
band-gaps UILj-lJ to both view points of their structural 
and dynamical properties. Crudely speaking, a nonlinear 
photonic band-gap excitation is an electromagnetic wave 
moving trapped to a defect in the dielectric material with 
an exponentially localized shape. 

Our aim in this work is to suggest a new approach to 
the global nonlinear feedback resulting from the super- 
position of the optical nonlinear properties of each sin- 
gle dielectric layer in a multilayer structure. In general, 
the propagation of electromagnetic field in a single di- 
electric layer is governed by the nonlinear Schrodinger 
equation(NLSE) i.e. 0: 



obtained as: 
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Solutions of this equation which take account of the local- 
ized feature of the light due to a local nonlinear feedback 
are stationary solitons. For a self- focusing nonlinear feed- 
back, the stationary soliton is bell-shaped(pulse) and is 



(z, t) = u(z) exp(—iu>t) , 
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u(z) = u a sech(\f~L) z) , u 2 D 
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In the spirit of the NLSE, this solution describes an 
harmonic time oscillation hidden in a spatial shell with 
a nonlinear topological shape. If we assume 10 to carry 
the optical information, then by virtue of its topological 
feature the nonlinear spatial excitation will act as a shell 
protector. To a mathematical viewpoint, (J2J is a spa- 
tial topological soliton with an internal frequency hence a 
breather soliton. Refering to the NLSE and from the ana- 
lytical expression of the soliton waveshape, to will depend 
on the intrinsic parameters of the breather. Thus @ is a 
gap soliton, where u is the amplitude of the soliton gap. 
If we now assume a superlattice resulting from peri- 
odic arrangement of single dielectric layers obeying equa- 
tion QJ, gap solitons of each of these single dielectric ma- 
terials become their foundamental modes. However, ow- 
ing to the broken translational invariance from the cou- 
pling of multiple layers, these foundamental modes will 
no more obey single NLSEs and the dynamics of a given 
n th optical component is from now on governed by the 
following equation: 
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We will assume all dielectric materials to be identical 
and denote their common foundamental frequency by 
VL . The key point of our study is the assumption that 
equation © is self-consistent, i.e. is valid for all optical 
components. Namely, the limit N — m should repro- 
duce the single-pulse soliton J2J. In this context, the 
sum in m runs over the entire superlattice and therefore 
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takes into account all foundamental modes. Doing so, 
the problem turns to that of a single dielectric layer in a 
potential bath(hereafter denoted "mean- field potential") 
created by its surroundings. This mean-field potential is 
adiabatic in the sense that it is switched on after soli- 
ton defects have formed on all the N dielectric layers. 
To allow for the interaction among solitons, we need to 
separate their centres of mass from their spatial coor- 
dinates on layers. In fact, this amounts to introduce a 
pinning-center in the overall coordinate of each founda- 
mental soliton so that the isolated single-pulse is readily 
rewritten: 
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With this more appropriate shape, equation J3J leads to 
the following equation for the spatial wavefunction of the 
light wave passing through a given n th optical component 



and interacting with an m one via a coupling J m - 
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Below, nonlinear solutions of this equation are discussed 
assuming two different contexts. We first consider the 
case of a uniform periodic arrangement of dielectric 
layers (regularly distributed single-layer feedbacks), and 
next the case of a dispersive coupling(dispersion-induced 
localized single-layer feedbacks) between optical compo- 
nents. 

For Jm-n — Jo where J is a constant, the mean-field po- 
tential experienced by the n th optical component is the 
sum over sech 2 functions. This sum is exact if we note 
the following relevant identity (£ = m — n): 
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where K, K' and E' are Complete Elliptic Integrals and 
sn is the Jacobi elliptic function of modulus v and pe- 
riod L. This relation follows from simple transformations 
of the usual Fourier series representation of Jacobi El- 
liptic functions [?J, taking into consideration the extra 
argument z n . The above identity, given in terms of the 
intensities of the local foundamental light waves, mim- 
ics a periodic lattice of single pulses mutually trapped 
at their centre-of-mass z n = nL where L is the sepa- 
ration between two interacting pulses. L increases as v 
is increased so that one can monitor the degree of com- 
pression of single pulses within the potential, namely by 
tuning the thickness d v of the mean-field potential wells 
which according to © is proportional to the width of the 
single-pulse defect(^=). On fig. ^ we sketch the space- 
dependent part of the global feedback potential(i.e. the 
sn 2 function). When v = 1 the mean-field potential re- 
duces to sech 2 \/TT ~z. In connection with this finding, 
below we will also see that we recover the original single- 
pulse soliton. Remark that in this limit the separation L 
between the single pulses is infinite. 
Replacing the identity © in equation (JSJ we find: 
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Equation Q, so-called first-order Lame equation 0, is 
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FIG. 1: The global nonlinear feedback potential. 



the actual one governing the spatial waveshape of the 
bunch of individual solitons passing through the n th op- 
tical component of the periodic lattice. This equation 



possesses three distinct bound states listed as 
ui(x) = u i(v)dn(x), 
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u 2 {x) = u o2 {v) sn{x), 

uji = n o K' [(1 + k 2 )K' - 2E'] f ~j . (9) 



FIG. 2: Spatio-temporal waveshape of the dn gap-soliton 
lattice. 



FIG. 3: Spatio-temporal waveshape of the sn gap-soliton 
lattice. 



u 3 (x) = u o3 (v) cn(x), 
uj 2 = tt K' [K' - 2E'\ 



(10) 



In general, these bound states form a complete orthonor- 
mal set and are normalized accordingly. However, to 
avoid a divergent light intensity we bind their amplitudes 
to the amplitude of the single-soliton gap u . The nor- 
malization condition then becomes: 
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and leads to: 



Ui(x) uj(x)dx = UgSij, L 1 = — 
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We can check that in the limit v — > 1 these bound states 
reduce to: 



u\(z) = u%{z) = u sech-\ 
U2(z) — 0, u>\ = 0. 
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On figs. |21 El and E] we display the spatio-temporal 
waveshapes^ = 0.98) of the three bound states. 

According to their asymptotic behaviours, the first and 
third bound states are two degenerate periodic gap soli- 
tons. The second bound state, which gives rise to a 
vanishing single kink in the limit L — > oo, is the state 
of broken translational symmetries of individual solitons 
due to the correlations of their (non-zero)centre-of-mass 
coordinates. The vanishing of this mode in the limit 
v — > 1 (where the separation between single pulses be- 
comes infinite) is a clear evidence of restored translational 
invariance. 

Let us now turn to the case of a dispersive coupling. 




FIG. 4: Spatio-temporal waveshape of the cn gap-soliton 
lattice. 



We will assume the simplest case of an exponentially de- 
creasing interaction between single dielectric layers i.e. 

Jm—n d Q C 

-A|tn-»| ^ 

Since the problem cannot be 
efficiently treated by anlytical means, we followed a nu- 
merical approach. Figs. El El 13 and El are shapes of 
the light intensity on a single layer within the multilayer 
structure, simulated assuming four different values of the 
separation L between N — 100 single pulses. The pulse 
width is 10 units correponding to a foundamental fre- 
quency fi ~ 0.01, while the soliton-lattice frequency 
Lu n = 10il o . A is set to 0.5 which we expect to ensure 
relatively long distance effects. The most stricking fea- 
ture in these figures is a less and less periodic distribution 
of single pulses in the gap-soliton lattice as their separa- 
tion L decreases. The gap-soliton lattice is clearly seen 
to transform into a bunched soliton state showing strong 
localization on single layers as solitons in the bunch get 
more and more compressed. Remark that fig. ?? corre- 
sponds to an inter soliton separation which is equal the 
single-pulse width. Actually, the behaviours emerging 
from the figures are not surprising since it is a natural 
fact that a decrease of L will enhance sensibility of the 
soliton shape to dispersive interlayer interactions. 
As a concluding remark, the new analytical gap-soliton 
derived above are artificial objects consisting of a bunch 
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FIG. 5: Interlayer-dispersion-induced localization of the gap- 
soliton lattice: L = 40. 



FIG. 8: The gap-soliton lattice for L — 5. 
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FIG. 6: The gap-soliton lattice for L = 20. 



of single localized gap-solitons which are compeled to 
share the same propagation medium, as opposed to the 
well-known periodons of finite single dielectric layers see 
e.g. Lidoriskis et al. 0. Therefore, they can be looked on 
as compressed single gap-solitons where the global mean- 
field feedback plays role of a compression potential. Since 
this global feedback forms from foundamental modes of 
each single dielectric layer, the thickness of its wells and 
amplitude of the potential barriers will be fixed by the 



single-soliton width and gap, respectively. That is to 
say, one can monitor the soliton compression within the 
global feedback. According to our analytical soliton so- 
lutions, the waveshape acquired by the gap soliton pass- 
ing through each dielectric layer reflects features of the 
compression potential, thus each of these layers actually 
convey a bunch of compressed single pulses. Within the 
bunch state the soliton gap is lowered, decreasing with a 
decrease of the modulus v. Recall that a decrease of v 
brings single solitons closer and closer. As the dispersive 
inter-layer interaction just sharpens the soliton-lattice on 
single layers, we also expect similar features (namely a de- 
crease of the soliton gap) to show up in this last context 
for a relatively large number of interacting single dielec- 
tric layers. 
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FIG. 7: The gap-soliton lattice for L = 10. 
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